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Abstract 

We present a measurement of the B°-B° mixing parameter Am^ using neutral B meson pairs 
in a 29.1 fb _1 data sample collected at the T(45) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. We exclusively reconstruct one neutral B meson in the semilep- 
tonic B° — ► D*'~i + v decay mode and identify the flavor of the accompanying B meson from its 
decay products. From the distribution of the time intervals between the two flavor-tagged B meson 
decay points, we obtain Am^ = (0.494 ± 0.012 ± 0.015) ps _1 , where the first error is statistical 
and the second error is systematic. 

PACS numbers: PACS numbers: 13.20.He, 11.30.Er, 12.15.Hh, 14.40.Nd 
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B°-B° mixing plays a unique role in the determination of basic parameters in the standard 
model (SM) of elementary particles. It is characterized by the oscillation frequency Am^, 
which is the difference between the two mass eigenvalues of neutral B meson states. In the 
SM, the mixing is due to second-order weak interactions known as box diagrams [l|] whose 
amplitudes involve V t d, sua element of the quark mixing matrix governing transitions between 
the top and down quarks [fj]. The mixing also induces large time-dependent CP violation 
in neutral B meson decays, which has been observed recently 0, ||. For such CP violation 
measurements, precise Am^ measurements are important. 

In this Letter, we report a Am^ measurement with 31.3 million B°B° pairs, collected 
with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) collider || 
operating at the T(AS) resonance. The time evolution is described as e~' Ai ' //T s°/(4rBo){l ± 
cos(Am^At)}, where the plus (minus) sign is taken when the flavor of one B meson is 
opposite to (the same as) the other, r#o is the lifetime of the neutral B meson and At 
is the proper time difference between the two B meson decays. At KEKB, the T(4S) is 
produced with a Lorentz boost of (3^ = 0.425 nearly along the electron beamline (z). Since 
the B mesons are approximately at rest in the T(4S f ) center-of-mass system (cms), At 
can be determined from the displacement in z between the two B decay vertices: At ~ 
(z rec - z taug )/(3jc = Az/p-fc. 

The Belle detector || is a large-solid-angle spectrometer that consists of a silicon vertex 
detector (SVD), a central drift chamber (CDC), an array of aerogel threshold Cerenkov coun- 
ters (ACC), time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located outside of the coil is instrumented 
to detect K° L mesons and to identify muons (KLM). 

We use the decay chain B° -> D*~i + u, D*~ -»■ D°n-, and D° -> K + 7r~, K + n-n° or 
K + 7t~7t + 7t~ [0. The large branching fractions and distinctive final states of the semileptonic 
decay allow for the efficient isolation of a high-purity B° sample. The event selection criteria 
are almost the same as those for our previous CP violation measurement ||. Charged 
particles are selected from tracks with associated SVD hits. Track momenta for D° — > 
K + n^7Y + n^ decays are required to be larger than 0.2 GeV/c. Candidate 7r° — > 77 decays 
are pairs of photons with energies greater than 0.08 GeV that have an invariant mass within 
0.011 GeV/c 2 of m n o and a total momentum greater than 0.2 GeV/c. We require the 
invariant mass within 0.013 GeV/c 2 of m D o for D° — ► K + tt~ or K + n~n + n~ , and —0.037 < 
M K + n - n o - m D o < +0.023 GeV/c 2 for JJ° — > K + ty'tt . The mass difference between the 
D*~ and D° candidates, M^, should be within 1 MeV/c 2 of the nominal value. The cms 
angle between the D*~ candidate and a lepton (an electron or a muon that has a cms 
momentum within 1.4 < p c e ms < 2.4 GeV/c) is required to be greater than 90 degrees. 
The energies and momenta of the B meson and the D*£ system in the cms should satisfy 
M 2 = (E c ™ - E c ™\) 2 - |p B cms | 2 - |p D Tl 2 + 2 IPB ms | \Pdh\ cos6 b ,d*£, where M v is the neutrino 
mass and 9b,d*i is the angle between pg ms and p^f. We calculate cos9b,d*£ setting M v = 0. 
Figure shows the cos 6 b,d*£ distribution. The signal region is defined as | cos 9 B,D*e\ < 11. 

We identify the flavor of the accompanying B meson from the properties of the decay 
products 0. Several categories of well measured tracks that have a charge correlated with 
the b flavor are selected: high momentum leptons from b — > c£~T7, lower momentum leptons 
from c — > s£ + v, charged kaons and A baryons from b — > c — ► s, high momentum pions 
originating from decays of the type B° — > D^*'~X (where X = 7r + ,p + , a^,etc), and slow 
pions from D*~ — > D°tt~ . Information extracted from each track is combined for the 6-flavor 
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FIG. 1: The cosOb^d'I distribution for the D*~t r v candidates. The circles with errors show the 
data. The solid line is the fit result. The total background and the D**£v component are shown 
by the dashed line and the hatched area, respectively. The inset shows the same figure with a 
logarithmic scale for the vertical axis. 

determination, taking into account correlations in case tracks in more than one category are 
present. For each flavor decision, we assign a MC-determined flavor-tagging dilution factor 
r, which ranges from r = for no flavor discrimination to r = 1 for unambiguous flavor 
assignment. It is used only to sort data into six intervals of r, according to estimated flavor 
purity. More than 99.5% of the events are assigned a non-zero value of r. 

We reconstruct the B° — > D*~£ + v decay vertex using the D° trajectory, the lepton track 
and the interaction-point profile (IP) convolved with the finite B flight length in the plane 
perpendicular to the z axis (21 /im). The reduced \ 2 °f the vertex is required to be less 
than 15. The method of reconstructing the tagging side B vertex is described elsewhere [|J. 

We find 16397 candidates after flavor tagging and vertex reconstruction. The signal 
fraction is estimated to be 80.4%. The backgrounds consist of fake D* mesons (7.8%), 
B — > D**iv events (7.4%), random combinations of D* mesons with leptons with no angular 
correlation (2.6%; called "uncorrected background") and continuum events (1.8%). Here 
D** consists of non-resonant D*ir components and charmed mesons heavier than D*. The 
background due to a combination of a fake lepton and a true D* from the same B meson 
is estimated with MC to be negligible. We estimate the fake D* background fraction from 
the D° mass sideband events and from fake D* events reconstructed with wrong-charge slow 
pions. The uncorrelated background fraction is evaluated by counting candidates where we 
invert the lepton momentum vector artificially. We estimate the continuum background frac- 
tion by scaling the off-resonance data (2.3 fb _1 ) with the integrated luminosity. We fit the 
cos9b,d*£ distribution in a range —10 < cos9b,d*£ < 11 to estimate the B — > D**iu back- 
ground fraction; the cos 9b, dh shapes for the signal and B — ► D**tv are modelled using MC 
and all the other background fractions and distributions are fixed from the aforementioned 
special background samples. 

The At resolution function for the signal, -R sig (At), is expressed as 



R sig (At) = g 1 G(At;^ 1 ,a 1 ) + (l-g 1 )G{At;^ 2 ,a 2 ) 
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where G(x; //i(2), 01(2)) is the main (tail) Gaussian component with / u 1 ( 2 ) and o"i( 2 ) as the 
mean and standard deviation, respectively, g\ is the fraction of the main component, Sip) 
is a scale factor that corrects our imperfect error estimation, and cr rec ( ta g) is the At error 
calculated from the vertex error for the reconstructed (tagged) B meson determined for 
each event. We extract the above parameters from the At distribution of the candidate 
B° — > D*~£ + i> events without distinguishing between different flavor assignments. The signal 
probability density function (PDF) is given by F sig (At) = A(At'; r B o)R sig (At- At')dAt', 
where A(At;r B o) = exp(— |At|/r B o)/(2r B o). We define the likelihood value for each event 
as Li = (1 — /bg)F S ig(Ati) + /bg-Fbg(Aij), where /b g is the overall background fraction of 
0.196 and Fb g (Atj) is the background PDF given by J2k fkF k (Ati). Here F k and f k are the 
PDF and the fraction, respectively, for each of the four background components. We use 
the signal PDF for the B — ► D**iv component. For the other background components, 
we use F fc (At) = /^[(l - f Sk )A(Af; r k ) + f Sk 5{At')]G{At - At'; n k , a k ))dAf, and a k = 



Sky c"rec + °fag- Here 5 (At') is the Dirac's delta function that accounts for components with 
small or zero lifetimes. The parameters in F k (At) are obtained from the upper sideband 
(0.155 < Mdifr < 0.165 GeV/c 2 ) for fake D*, the lepton-momentum-inverted events for 
uncorrelated background, and off-resonance data for continuum background, respectively. 
We perform a likelihood fit to determine R sig (At) with t b o = 1.548 ps || and with the 
background parameters fixed to the obtained values. We find that the fraction of the main 
component is large (gi = 0.87ig;og) and the At error estimation is correct (Si = 0.99to;?o) • 
The typical rms resolution on At is 1.43 ps. 

We determine Airtd by an unbinned maximum-likelihood fit to the At distributions. We 
define the likelihood value for each event as follows: 

4 OF(SF) = (1-4){(1-W)C (SF) (A^) 
+ f D ** e „FSif\AU)} 

, f l ST f l fOF(SF) p OF(SF), A , x 
+ Jhgl^JkJlk b k ( A *iJ) 
k 

where OF (SF) denotes B°B° (B°B° or B°B°), i.e. a state with the opposite (same) flavor, 
F°f (SF) , F%lfP and F° F(SF) are PDFs for the OF (SF) signal events, D**£v decays and 
other backgrounds, respectively, (I = 1, 6) is an overall background fraction excluding 
B — ► D**tv in each r region and (1 — fl g )fD**iv corresponds to a fraction of D**tv decays. 

Other background fractions fi^ F ^ SF \ where the relation + /;| F = 1 holds, are obtained 
from the control samples for the uncorrelated and fake D* backgrounds, and from MC events 
for the continuum. 

The PDFs for the OF and SF signal events are given by 

F°f (SF) (At) = [™ V° F f F) (At')R sig (At-At')dAt>, 

J — oo 

where 

vZ F F F) (At) = — {1 ± (1 - 2wi) cos(Am rf At)}. 

At b o 

Here wrong tag fractions wi are also determined simultaneously. The PDF for B — ► D**iv 
is given by a sum of B° and B + components as F®f*f F '(At) = (1 — f B +)F® g ^ SF \At) + 
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f B +F B + ^ SF ^(Ai), where f B + is the B + fraction in the B — > D**£u background. The B + — > 
D**lv background PDFs F B + (SF) (At) are given by the following functions convolved with 
# sig : V°l = (l-w l B+ )V B + andVf + = w l B+ V B+ , where TV (At) = {l-rsf s )A{At]r B +T , B+ ) + 
rgfsS(At) and w l B+ is the wrong tag fraction determined from a B + — > D°n + sample. Since 
the fake D* background includes a mixing component, we use a function that has the same 
form as the PDF for B — > D**lv with parameters obtained from the M^m sideband events. 
Other background PDFs do not distinguish between SF and OF events and are the same as 
those used to determine the resolution function. 

We perform the fit with 10 free parameters (listed in Table |) to the At distributions 
of SF and OF events in the signal region and the B — > D**£u dominant region defined as 
— 10 < cos9b,d*i < —1.1. In this way, the background parameters f B +, fs and t' b+ are 
determined simultaneously. Additional correction factors, r$ and r B +, are introduced only 
in the signal region to account for the difference between the two regions. We use MC to 
determine r s = and r B + = 1.04 ± 0.02. 

The fit result is summarized in Table |. Figure [| shows the observed At distributions 
for the OF and SF events. Figure |3| shows the corresponding flavor asymmetry, A(At) = 
[N OF (At) - N SF (At)]/[No F (At) + N SF (At)], where N OF (sf) denotes the number of OF (SF) 
events. 

TABLE I: Summary of mixing fit. Errors are statistical only. For each wrong tag fraction, the r 
interval and the number of candidate events are also shown. 

parameter result 
Am d 0.494 ±0.012 ps" 1 

wi (0 < r < 0.25, 6360 events) 0.467 ± 0.010 

w 2 (0.25 < r < 0.5, 2364 events) 0.360 ± 0.016 

w 3 (0.5 < r < 0.625, 1453 events) 0.254 ± 0.020 

w 4 (0.625 < r < 0.75, 1702 events) 0.182 ± 0.017 

w 5 (0.75 < r < 0.875, 1958 events) 0.103 ± 0.014 



w 6 (0.875 < r < 1, 2560 events) 0.032 ± 0.010 

hO.12 
-0.13 



f B + 0.70+ ' 12 



fs 0.28±g:Sg 
r' B+ 1.87lg:SS P« 



The systematic errors are summarized in Table [TT]. D** branching fractions used in this 
analysis are based on theoretical assumptions 0. We set each such branching fraction to 
unity in the MC (with all others set to zero), and repeat the analysis; we take the largest 
variation on the Am^ result as the systematic error. To account for uncertainties in the tails 
of the vertex resolution, we measure Amd by setting the upper limit on | At\ that ranges from 
5 to 55 ps. We take the largest difference from the main result, which is obtained without 
the | At | upper limit, as a systematic error. Systematic errors from the background PDFs 
are obtained by varying each shape parameter individually, repeating the fit procedure, 
and adding each contribution in quadrature. We also perform a MC study where we obtain 
background PDFs by two methods: one from the background control samples, and the other 
directly from the signal region. The difference between two Amd fit results is included in the 
systematic error. A fit with the B — > D*lv and B — > D**£u MC events yields the Amd value 
that is consistent with the input value within 1.7a. We conservatively take the difference, 
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FIG. 2: The At distributions for (a) the OF events and (b) the SF events. The solid lines are the 
result of the unbinned maximum likelihood fit. The dashed lines show the background distribution. 
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FIG. 3: The observed time-dependent flavor asymmetry. The curve is the result of the unbinned 
maximum likelihood fit. 



which we attribute to MC statistics, as a systematic uncertainty. The systematic error due 
to the IP constraint is estimated by varying (±10/xm) the smearing used to account for the 
B flight length. Other sources of systematic errors are obtained by changing each parameter 
by la, repeating the fit procedure and adding each contribution in quadrature. We also 
perform a Am^ fit in each r region and find no systematic trend. 

In summary, we have measured the B°-B° mixing parameter Am,; using B° — > D*~£ + v 
decays in a 29.1 fb" 1 data sample collected with the Belle detector at the KEKB e + e~ 
collider operating at the T(4S f ) resonance. From an unbinned maximum likelihood fit to the 
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TABLE II: Summary of the systematic errors (ps 1 ) on the Amj measurement. 



source 


error (ps 1 ) 


D** branching fractions 


0.007 


At range 


0.007 


Background shape 


0.006 


Resolution function 


0.006 


B° lifetime 


0.005 


Fit bias 


0.004 


Background fraction 


0.003 


B -> D**£v fraction 


0.002 


IP constraint 


0.002 


wrong tag fraction 


< 0.001 


B^ shape parameter 


< 0.001 


total 


0.015 



At distributions for B pairs with the same and opposite flavors, we obtain 

Am d = 0.494 ± 0.012(stat) ± 0.015(syst) ps -1 . 

The result is one of the most precise measurements performed so far, and is consistent 
with the world average value of Am^ = 0.472 ± 0.017 (ps -1 ) M as well as other recent 
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